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Abstract 



The nature of electrostatic interactions involving polyanions modulate the proper- 
ties of both synthetic and biological macromolecules. Although intensely studied for 
decades the interplay of many length scales has prevented a complete resolution of 
some of the basic questions such as how the electrostatic persistence length (Ig) varies 
with ionic strength (/). In this review we describe certain characteristics of poly- 
electrolytes (PAs) and polyampholytes (PAs), which are polymers whose monomers 
have a random distribution of opposite charges. After reviewing the current theoret- 
ical understanding of the dependence of Ig on I we present experimental data that 
conform to two distinct behavior. For RNA and DNA it is found that that Ig ~ 
whereas for some proteins and other polyelectrolytes Ig ^ A scaling type theory, 
which dehneates charge correlation and pure polyelectrolyte effects for the shape of 
PAs that is valid over a wide range of salt concentration is described. We also use 
theory and simulations to argue that the distinct stages in the kinetics of collapse of 
PAs (with a net charge that is small enough to induce globule formation) and PEs 
(relevant for RNA folding) are similar. In both cases the major initial conformation 
change involves formation of metastable pearl-necklace structures. In the coarsen- 
ing process large clusters (pearls) grow at the expense of smaller ones by a process 
that is reminiscent of Lifshitz-Slyozov mechanism. Finally, recent theories and single 
molecule experiments on stretching of single stranded DNA and PEs further sheds 
insights into the complex behavior of charged macromolecules. The survey, which is 
limited to very few topics, shows the importance of polyelectrolyte effects in a wide 
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4.1 Introduction 



The charged nature of biological macromolecules (not limited to DNA, RNA, and 
proteins) and synthetic polymers makes the study of polyelectrolyte (PE) properties 
important. Shape fluctuations in charged biomolecules control a number of func- 
tions such as transcription, packaging of DNA and RNA into phage heads as well 
as interactions involving nucleic acids with other biological macromolecules (higher 
order structural organization of DNA upon interaction with histones for example). 
An example of functional significance is the spectacular phenomenon of reversible 
condensation of DNA into toroidal structure in the presence of multivalent cations. 
Upon condensation the volume occupied by DNA decreases by nearly four or more 
orders of magnitude. Description of counterion-induced condensation of DNA and 
synthetic PEs requires accounting for density fluctuations associated with counteri- 
ons and their coupling to conformations of the macromolecule. Similarly, the phase 
diagram of charged synthetic polymers, especially in aqueous solution, is complicated 
and is critically dependent on ion concentration as well their size and shapes. The 
behavior of PEs is governed by a number of factors such as the intrinsic properties 
of the backbone (extent of flexibility), length scale associated with ion fluctuations, 
and even the chemical nature of the counterions. As a result it has been difficult to 
develop a coherent theory that can capture all the properties of polyelectrolytes. 

In this Chapter we confine ourselves predominantly to the study of properties of an 
isolated PE and not so unrelated problem of polyampholytes (PAs) and develop con- 
cepts that have been particularly useful in understanding a number of aspects of RNA 
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folding and a few problems associated with DNA. One focus, described in Section 2, 
is to elucidate the current understanding of how the persistence of charged chains, 
a basic property of polymers, depends on salt concentration. Although considerable 
work has been done in this area (see an excellent reviewl^ by Barrat and Joanny that 
still remains relevant) a clean solution to this problem has not been obtained. This 
is probably because there are multiple scenarios that depend on the molecular weight 
of the polymer and the concentration of chemical characteristics of the counterions. 
In Sections 3 and 4 we provide a simple picture of shapes of polyampholytes (PAs), 
which are polymers with random distribution of positively and negatively charged 
monomers. In Section 5 we describe the response of PEs to mechanical force, which 
have provided insights not only into the problem of electrostatic contribution to per- 
sistence length of charged polymers but also have been instrumental in describing 
length scale-dependent elasticity of single and double stranded DNA. We describe 
simulations (Section 6) and theory (Section 7) of collapse of PAs and PEs, which 
prove to be important in understanding folding of RNA. In Section 5 we describe the 
response of PEs to mechanical force, which have provided insights not only into the 
problem of electrostatic contribution to persistence length of charged polymers but 
also have been instrumental in describing length scale-dependent elasticity of single 
and double stranded DNA. The Chapter ends with concluding remarks in Section 8. 
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The problem of electrostatic persistence length of polyelectrolytes has been studied 
for more than three decades using a variety of methods, including analytical treat- 
ment, experiments and simulations. Due to an interplay of a number of length scales, 
compared to the case of neutral polymers, a complete theoretical understanding of 
persistence length changes in PE as a function of salt-concentration is difficult. The 
electrostatic persistence length is defined as the contribution to the chain total per- 
sistence length in addition to the "bare" (intrinsic) value by electrostatic interactions 
among the charged monomers. Theoretical description of the characteristics of an iso- 
lated PE chain began with two pioneering and independent studies by Odijk and by 
Skolnick and Fixman (hereafter referred to as OSF theory), which addressed the prob- 
lem of electrostatic contribution to persistence length of an intrinsically stiff chain. 
Besides the usual chain intrinsic stiffness of a neutral chain, which is characterized 
by bending rigidity (persistence length, Kuhn length, monomer length, interaction 
strength), electrostatic interactions among monomers can further stiffen the chain. 
The interplay of the chain intrinsic stiffness and the electrostatic interactions, which 
are usually modeled by screened Coulomb potential using Debye-Hiickel theory, po- 
tentially gives rise to a number of scenarios, which require a number of approaches. 
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4.2.1 Intrinsically stiff chains: The Odijk-Skolnick-Fixman 
theory 

Odijk'2' and, independently, Skolnick and Fixmanl^ were the first to develop a theory 
for persistence length of polyelectrolyte. Following Barrat and Joanny'I' we model a 
polyelectrolyte as a worm-like chain with a contour length L and "bare" persistence 
length Iq, and that carries charges separated by a distance A along its contour. The 
electrostatic interaction between two charges separated at a spatial distance r is 
modeled by the Debye-Hiickel potential: 



exp ( 

VDH{r) = ksTlB (4.1) 



where Ib = e^/{ekBT) (in Gaussian units as in Refs.^'^) is the Bjerrum length and 
is about 0.7nm at room temperature, and Ad (or, interchangeably, is the Debye 
screening length, which is given by 

with / being the ionic strength of the solution. 

For simplicity we confine the chain in two dimensions so it is a planar curve. The 
results in three dimensions are similar. The chain configuration is characterized by 
the angle function (of arc-length s) 9{s) = t(s)t(0) where —L/2 < s < L/2 and t{s) 
is the tangent vector at s. We assume that the chain is almost linear (strong charge 
limit) or 6{s) is small, so that the electrostatic interaction between two charges at 
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si and s2 can be expanded around the rod-like conformation, VDHi\s2 — si\). With 
these conditions the chain Hamiltonian H{6) can be written for a given configuration 

{e{s)} as: 



H[e] =Ho + hsT ds ds''^ [m{s - s') + K{s, s')] "^^^''^ 

2 J~L/2 J-L/2 ds 



L/2 J-L/2 ds 



(4.3) 



Here, Hq is the electrostatic energy of a rod, the term proportional to Zq is the usual 
bending energy of the worm-like chain that is proportional to the square of radius of 
curvature, and the electrostatic interactions described using the kernel K{s, s') which, 
for s > s', is 



K{s, s') = j-f^ ds, ds2v'^^{s2 - s,) ^'\ ''^ (4.4) 

J-L/2 Js [S2 - Si) 



with v'j-)^{s) = The result has been obtained with the only assumption that, 

within Xd the chain remains in an almost rod-like configuration. For very long chains. 



L/Xd ^ 1, we can extend the limits in the integrals of (4.4) to infinity and K[s, s') 
becomes a functions of s — s', 



1 f°° 
bA^ Jq X + s 



(4.5) 



which has the Fourier transform 



oo 



K{q) = / ds exp{iqs)K{s) 
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where 

losF = Ib/{^A\^) (4.7) 

is called the celebrated Odijk-Skolnick-Fixman (OSF) length. The statistical proper- 
ties of the chain are obtained by integrating the Boltzmann factor exp {H{[6]/kBT) 
over all possible configurations, i.e. over all functions 0{s). Since the energy is a 
quadratic function of 6, the integration can be carried out analytically. For instance, 
we can calculate the mean squared angle (6'(s)^) between the chain directions at the 
origin and at arc-length s. This quantity characterizes the local flexibility of the 



chain. With the approximate Hamiltonian (4.3), we have 



4 sm\qs/2) 



TT 



im') = - dq 7—^. (4.8) 



q' /O + K{q) 



For a neutral chain, K{q) = in Eq. (4.8) and we have 



{0{sf) = (4.9) 



which varies linearly with s, and the proportional constant is the inverse of the per- 
sistence length. 

For a charged chain, (^(s)^) cannot be analytically calculated. Thus, to describe 
the expected scaling behavior, we work in the limits of small and large s. At large s 



(4,10) 
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which expresses the fact that at long length scales, the chain conformation can be 
described by an effective persistence length IO + Iqsf, which is the sum of a "bare" and 
of an electrostatic contribution. It indicates that the influence of the screened elec- 
trostatic interactions can extend far beyond their range A^, since in weakly screened 
solutions, losF is much larger than \d. The electrostatic persistence length also de- 
creases with salt concentration as Iqsf — I^^-, whereas the Debye screening length 
has a slower decay \d — /^^^^ 



On the other hand, for small s, (4.8) reduces exactly to that of a neutral chain in 



Eq.(4.9). This means that the chain statistics at short scales are not modified by elec- 



trostatic interactions. The crossover between the "intrinsic" regime described by (4.9) 



and the "electrostatic" regime described (4.10) takes place when the electrostatic in 



teractions become strong enough to perturb the statistics of the neutral flexible chain. 
The crossover length Sc can be obtained qualitatively from the following argument. If 
a small chain section, of length s < k^^, is bent to form an angle 6, the cost in "bare" 
curvature energy is ksTlOO^ /s, while the electrostatic energy is kBTlB^s / A^iO"^ / s). 
The two energies are comparable for s ~ Sc, which gives Sc ~ A(/0//b)^^^. 

The picture that emerges from this calculation is that the chain flexibility depends 
on the length scale. At short scales, s < Sc, the chain structure is determined by its 



bare rigidity 10, while the electrostatic rigidity (4.7) dominates at large scales. Large 
length scale properties, such as the radius of gyration can be determined by applying 
the standard formula for semi-flexible chains of persistence length Iqsf- 

The only approximation required to obtain the Odijk-Skolnick-Fixman length is 



the expansion that yields equation (4.3). The calculation is therefore consistent if the 
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angle {6{k,~^)'^) is small compared to unity, which holds good if the bare persistence 
length is much larger than Iqsf- In other words, the requirement {9{k^^Y) ^ 1 is 
then equivalent to ScIqsf / + Iqsf)) — Sc/lO ^ 1. i.e. the angular deflection 
of the chain must be small when the crossover region is reached. In other words, 
the angular fluctuations that take place before the electrostatic interactions can come 
into play and rigidify the chain should not be too strong. If these fluctuations are too 
large, i.e. if the chain is too flexible, the perturbation expansion that underlies the 
OSF calculation can break down. The criterion for the validity of the calculation can 
be simply written, in the limit of weak screening, as 10 > A^/Ib- This implies that the 
OSF calculation should be directly applicable to stiff chains such as DNA {10 ~ 50nm), 
but has to be reconsidered for flexible chains such as polystyrene sulfonate {10 ~ Inm). 

Multiple lines of inquiry have questioned the validity of the OSF theory (see 
below). All of them hinge on the assumption that when chain fluctuations are taken 
into account the persistence length is renormalized and the dependence on deviates 
from the predictions of OSF theory. Starting with the work of Barrat and Joanny,!^ 
several variational calculations have shown that, when the chains are intrinsically 
flexible, the electrostatic persistence length Zg ~ A^. This dependence is valid when 
the controlling parameter ~ 0(1), i.e., when the chain is in the so-called non- 
asymptotic regime. Li and Wittenf^ argued that, unlike the variational theories, 
approximate inclusion of fluctuations still leads to the result predicted by OSF even 
if the backbone is flexible. In other words, the OSF result is always valid at least for 
long chains regardless of the backbone stiffness. In the above referenced variational 
theories it is generally believed that the electrostatic interactions generally stiffen the 
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chain significantly so that the trial Hamiltonian should consist of terms that account 
for the electrostatic-mediated interaction rigidity. 

The controversies associated with the validity of the OSF theory were further 
addressed using two related theoretical studies that used different variational meth- 
Q^gJUSl Using the most general variational approach these studies showed that the 
OSF result is recovered whenever <^ 1. Thus, in both the stiff and intrinsically 
flexible limit the classical result is obtained provided the chain is long. Indeed, these 
theories supported a scaling-type argument due to Khokhlov and Khachaturianf^ 
(KK), who argued that chain fluctuations merely renormalize A by a blob with di- 
mension D with charge. The physical that emerges is that fluctuations do not 
invalidate the qualitative aspect of the OSF result. Rather, they effectively reduce 
the direct distances between charges so as to renormalize the bare parameters lo, A, 
charge per blob, as was first recognized by KK. In this limit, local fluctuations inside 
the blobs are strong enough to reduce significantly the direct distance between two 
consecutive charges by the factor which is much less than one. This results 

in much stronger effective Coulomb repulsion between two consecutive blobs. Since 
the length scale of Ip is much larger than the blob size D and ordering as implied in 
these cases does not refer to the local structure of blobs, we can expect that the local 
fluctuations inside a blob not to affect significantly the property of the chain much 
beyond the length scale D. As a result, we obtain qualitatively similar dependence 
of /p on Ad as that obtained for stiff polyions. Note that this is relevant only for the 
asymptotic case where each blob contains a large number of segments which roughly 
obey Gaussian statistics. The variational theories also predict interesting crossover 
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behavior in which the electrostatic persistence length ~ where the exponent y < 1. 

4.2.2 The linear and sub-linear dependence of /e on Ad 

The OSF theory has formed the basis for interpreting both experiments and simula- 
tions. Indeed, a number of experiments have confirmed the quadratic dependence of 
the electrostatic persistence length on the Debye length jSliUllS] However, many 
other experimentsl^^^ and simulations'^^'^ have shown that the dependence of le on 
Xd is linear or even sub-linear. The reason for this discrepancy seems to be because 
the OSF theory was originally derived for intrinsically stiff polyelectrolytes, those 
that have the intrinsic persistence length much larger than separation of charges on 
the backbone Iqsf ^ such as dsDNA and carboxymethylcellulose (CMC).'^In the 
latter experiments, intrinsically more flexible polyelectrolytes, such as ssDNA and 
RNA were examined. Thus, chain fluctuations or finite size effects can indeed pro- 
duce deviations from OSF theory. The deviations from OSF behavior have lead to 
other approaches to study flexible polyelectrolytes. 

Barrat and Joanny, on the other hand, used a variational method in which they 
replaced the actual Hamiltonian with a "trial" one Ht, and minimized the resulting 
free energy F^ar = (Ht) —TSt, where the average is taking with the Boltzmann weight 
exp [—PHf] with /3 being the inverse of the thermal energy ksT. Here Ht is that of 
a neutral chain under uniform tension. For flexible chains, they also found that the 
total persistence length scales as Xd. Dobrynin'231 also modeled a polyelectrolyte as a 
wormlike chain with electrostatic interactions and evaluated the bending angle fluctu- 
ations in the frame work of the Gaussian variational principle. For semiflexible chain 
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and strongly charged flexible chain, Ig ~ Xd, whereas the dependence is sublinear for 
weakly charged chains. Ha and Thirumalai (HT)'251 used a self-consistent variational 
theory to calculate the total persistence length of a polyion. The theory is general 
for both flexible and stiff chains. They found that 

L + ^OSF, if ^OSF ^ lo 

lp-< ^ ^ (4.11) 

(loUJc)^^^ Xd, if Iosf'^L 

Therefore, for stiff chains (/o ^ ^osf) the electrostatic persistence length le is Iqsf ~ 
Xjj, whereas for flexible chains {Iqsf ^ D, h is not much different from Ip and scales 
as X]j. 

The effect of valence in driving collapse of RNA was vividly illustrated using 
small angle X-ray scattering (SAXS) experiments on a ~200 nucleotide Azoarcus 
ribozyme.'^ From the measured scattering intensity they obtained the distance dis- 
tribution, P(r). The dependence of the square of the radius of gyration, which cor- 
responds to the second moment of P(r), depends sensitively on the valence of the 



counterion (see Fig. 4.1). The ribozyme is extended at low cation concentrations and 



is compact at elevated values of the counterion concentration [Fig. 4.1| . The collapse 



transition is highly cooperative in Mg^"*" and is much less so in Na"*" [Fig. 4.1 . It was 



found that P{r) (Fig. 4.2) could be well fit using the asymptotic form for the end-to- 
end distribution for WLC chains.'^ The experiments showed the persistence length of 
this RNA changes dramatically from about 3nm at low ionic concentration to about 
Inm at high salt concentrations. More importantly. For both Na+ and Mg^"*", the 
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persistence length scales as Xj^, which is consistent with the OSF theory. From this 
finding, we find that the intrinsic persistence length of RNA is ~ Inm. 



The evidence of the linear dependency of the electrostatic persistence length on 
A/) is provided in a number of experimental studies. Tricot^ analyzed the intrinsic 
viscosity-molecular weight dependence of a number of polyelectrolytes such as car- 
boxymethylcellulose in solutions with / ranging from 0.005M to l.OM. The linear 
dependency is found for other types of charged polymers, whose total persistence 



length is on the order of lOnm or less (see Fig. 4.4). Tinland et al.^ measured the 
self-diffusion coefficient of ssDNA fragments using florescence recovery after photo- 
bleaching to infer the persistence length. The fltted total persistence length measured 
in Angstroms is found to be Ip = 6.4210"^ + 4/"^/^ with I measured in molars. Per- 
haps, the most convincing evidence for the linear variation is found in recent single 
molecule stretching experiments. Saleh et al. combined single molecule stretching 
data of ssDNA at very low forces with scaling arguments to obtained Zg ~ /-o-5i±o.o4 
or If, ~ j-o.4o±o.o4 rjj^g similar behavior is once again reproduced using a general 
property of the force extension curve at point where the relative extension is about 
1/2 by Toan et al.'^ Toan et al.'^ also found the linear dependency for synthetic 
RNA by fltting the force extension curve data by Seol et al.'^ using the Thick chain 
model.™ 
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4.3 Polyampholytes 

In this section, we present a simple physical model for describing a polyampholyte 
(PA) chain, especially the weakly-charged case, with a random distribution of oppo- 
site charges, in an electrolyte solution (e.g., NaCl dissolved in water). Here, we do 
not attempt to elaborate on the known results in the literature.'^ There has been 
much progress in understanding the physical properties of a PA chain under various 
conditions. For instance, the shape of a PA chain has been studied extensively for 
varying solvent quality (good vs poor solvent) and for different charge distributions 
(e.g., a varying degree of excess charge) (see Refs.l^ and references therein). How- 
ever, the earlier effort has been focused on the low-salt and high-salt limits. Our 
main motivation here is to offer a unified picture, in which both limits are integrated 
coherently. In appropriate limits, our description reproduces known results. Also, in 
our approach, the crossover between the two limits is captured in a more transparent 
manner. This effort will be beneficial for further illustrating various competing effects 
in shaping a PA chain and for offering a more coherent picture of such a system. 

4.3.1 Charge correlations and screening 

Much of our discussion will rely on the concept of "screening," which is well-understood 
for a simple electrolyte (e.g., NaCl dissolved in water) .'^ Select any charge and place 
it at the origin. Since the entire system is electrically neutral, the charge at the 
origin will be surrounded and shielded by the ionic cloud of the opposite charge. If 
the energy of this system tends to shrink the ionic cloud,'^ the entropy opposes this 
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tendency. This means that the thickness of the ionic cloud depends on the ion concen- 
tration and the temperature.'^ This idea can be extended to PA charges. While both 
species (i.e., PA charges and salt ions) contribute to screening, their roles in shaping 
the PA chain are opposite. Along this line, it proves useful to distinguish between 
"self-screening" (screening of a PA charge by other PA charges) and salt-screening 
(i.e., screening of PA charges by salts). Here self-screening refers to the tendency of 
opposite charges on the chain to be spatially correlated. This effect is responsible 
for the electrostatic compaction of a PA chain (under the right conditions) and thus 
competes with the latter effect, since the surrounding salt ions tend to diminish this 
tendency;!^ there is a competition for the spatial correlation of a PA charge with 
salt ions. Despite the difference between self-screening and salt-screening, the elec- 
trostatic free energy of a PA chain can be obtained by treating both kinds of charges 
as forming a Debye-Hiickel solution (or an electrolyte at the low electrostatic coupling 
limit'^. This is most evident within the theoretical framework known as the random 
phase approximation.'^^'^ The level of approximations assumed for this is similar to 
that for the Debye-Hiickel approach. At this level, chain connectivity does not influ- 
ence the electrostatic free energy of such a (randomly-charged) PA chain (for a given 
chain size) .121122] Obviously, non-electrostatic terms are different for the PA chain and 
for the surrounding electrolyte. 

Consider a PA chain consisting of monomers of size b each, which is inside an 
imaginary volume V of radius R. The PA chain is assumed to be weakly charged. 
Let f± be the fraction of positively and negatively charged groups. Then /± ^ 1. 
For a neutral chain, /+ = /_ (= /). Let n± be the total density of positive and 
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negative ions whether they are free or on the chain; let e be dielectric constant of 
the solvent and ksT the thermal energy. The Bjerrum length, at which the Coulomb 
energy of two electronic charges is equal to the thermal energy kBT, is then expressed 
as Ib = /eksT. (Here we adopt Gaussian units as in Refs.'^'^ Finally, the Debye 
screening length, is given by the relation = ATrlB^n^ + Outside V, n± 
tends to a constant, hq (2no is the salt concentration at infinity). Inside V, however, 
the salt ion concentration can be perturbed by the PA charges. Nevertheless, we 
write as = /tp^ + Kq inside V, where ~ NJIb/R^ and = SirlBno- Our 
DH approach will reproduce known results based on a more elaborated approach, in 
which the salt concentration inside V is determined more systematically .^1 



The electric potential due to a point charge e is given by 



e e 



e r 



(4.12) 



This is a potential screened by both PA charges and salt ions. To isolate the con- 
tribution of PA charges to screening, one has to subtract the salt contribution as 
follows 



PA 



lim 

r-5>0 



e e 
e 



-Kor 



Kq). 



(4.13) 



The electrostatic energy gain per charge due to self-screening or the polarization free 
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energy of the PA chain is 



EpA 



2 ksT 




(4.14) 



The factor 1/2 is to avoid double counting of interaction pairs. The total polarization 
free energy of the PA chain can be written as 



where kpa is expressed explicitly in terms of PA parameters. Note here that the 
electrostatic free energy differs from the energy contribution by a numerical prefactor. 
Clearly, this free energy favors chain collapse. 

4.3.2 Flory theory: a G-solvent case 

In Flory theory, the non-electrostatic free energy can be expanded in powers of the 
monomer density p = N/ as 



where v and w are the second and third virial coefficient, respectively.'^ The total 
free energy is the sum of FpA and F: 




(4.15) 



F 



(4.16) 




(4.17) 
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Here we mainly focus on a ^ solvent in which v = 0, even though our general approach 
can readily be extended to other cases (see Refs-S^^. For f = 0, we show our approach 
reduces to existing ones in low salt and high salt limits. In this case, the two-body 



term, i.e., the first term in Eq. 4.16, vanishes. The equilibrium size of a PA chain 
is then determined by the balance between FpA and the three body term, i.e., the 
second term in Eq. 4.16 ^ The equilibrium R satisfies 



NfhV (Nfl 



V R^ 



B , 2 
+ /€0 



-1/2 



N 
R^ 



0. 



(4.18) 



Because of the functional form of this relation, one can easily find that R ~ N^^^. 
The dependence of R on other parameters is more involved. After analyzing the 



relation in Eq. |4.18| in the two limiting cases, — )■ and kq — )■ oo, we will present 
our numerical solution for the intermediate range of kq. For simplicity, we set w = 1. 



Low-salt limit: k — )■ 0. In this case, Eq. 4.18 reproduces the known resulfl^ 



R ~ N^/^b 



fl 



B 



1/3 



(4.19) 



The chain size decreases as / increases, as expected based on the following picture. 
For larger /, the effect of self-screening or the PA effect is more pronounced and 
shrinks the chain size more effectively. 



High-salt limit: k — )■ oo. In this case, Eq. 4.18 reduces to 



(Nfl 

[ R' J ^0 



i?3 



0. 



(4.20) 
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This reproduces the known result for the high-salt limit'^ 

.^.V3,(-J)'. ,,,, 

The chain is swollen by added salts. This finding is consistent with our earlier asser- 
tion that self-screening competes with salt-screening. 

Crossover. At the crossover between the two limits, the PA size expressed in 



Eq. 4.19 becomes comparable to that in Eq. 4.21 This leads to the following crossover 
conditional 

flB 



kq. (4.22) 



62 

In fact, the term on the left hand side is the reciprocal of the self-screening length, 
KpA, for a PA in the low-salt limit: 



'^PA~V^^~-^- (4-23) 



The crossover condition thus reads kpa ~ kq. 

For the intermediate salt range, we have solved Eq. |4. 18| numerically and plotted 



our results in Fig. 4.6 Our results in Fig. 4.6 show how added salts expand the PA 



chain. A key to understanding this behavior is that the tendency for PA charges to 
be correlated is weakened by the surrounding salt ions, since they are equally corre- 
lated with the salt ions. This competition underlies the Kq dependence of R shown 
in Fig. [46 
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In summary, our DH approach combined with Flory theory offers a unified frame- 
work for describing the equihbrium properties of a polyampholyte solution in the 
presence or absence of added sahs. It not only reproduces known resuhs in the low- 
salt and high-salt limits but also shows the crossover between the two limits. In the 
next section, we complement this approach using the blob picture of a PA chain. 

4.3.3 Blob picture of a polyampholyte chain 

The notion of blobs has been proven to be useful for illustrating various aspects 
of polymer systems.'^ The section is devoted to presenting a blob picture of a PA 
solution, offering a "visual" guide to the approach presented in earlier subsections. 
For a weakly charged PA chain, which we consider here, the chain statistics at short 
length scales will not be perturbed by the PA effect. Within this length scale, the 
PA chain resembles an ideal chain or a random walk (in a 6 solvent). Beyond this, 
the PA effect will collapse the chain. This means that there exists a special length, 
denoted as ^pa, that separates these two regimes. To obtain this length, consider a 
section of the chain, consisting of g monomers. For sufficiently small g {> 1), the size 
of this section is given by r ~ bg^^"^. The free energy of this section is dominated by 
the entropic term r'^/gb'^. As g increases, the PA term becomes more important. To 
set up the crossover condition, find g at which the two terms are equally important: 




(4.24) 
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This leads to 

9flB (\/^ + -o--o)-l. (4.25) 



The g value that satisfies Eq. 4.25 is the minimum number of steps the random walk 
has to take until it starts to feel the PA effect. The blob size is then ^pa ~ bg^^'^. For 
N ^ g, the chain can be viewed as a compact stack of many such blobs. This means 
that R ~ ^pa(A^/^)^/^ ~ N^/%g^/\ 



Low-salt limit. As kq — )■ 0, Eq. 4.25 results in (7 ~ b"^ /U'IbY and ^pa ~ b"^/ fls 
Kp\, consistent with ReiJ^ 



High-salt limit. For large kq, Eq. 4.25 reduces to (7 ~ t^yf' / [flsY ■ The length 



equivalent is then ^pa ~ t^ob'^ / U^bY ~ '^o/'^pa' agreement with Ref.^l 

In both limits, this blob picture reproduces the results for R in Eqs. 4.19 and |4.2l] 



obtained by free energy minimization. Also one can show that FpA ~ ksT x (N/g). 
This is simply the free energy cost for "redirecting" the random walk so as to fill the 
space compactly. Hence ksT per blob. 

Intermediate salt concentration. For the intermediate range of kq, g has been 



obtained numerically and plotted in Fig. 4.6 (see the dotted lines) for a few choices of 



/; we have used the same color scheme as for the R-kq plot. As shown in the figure, 
^PA grows in magnitude as kq increases. This is paralleled with the earlier finding 
that the PA chain is swollen by salt. 

Surface tension. Our DH approach to a PA chain is based on the assumption that 
N is arbitrarily large. In this case, the "surface effect" is minimal, as long as the 
PA chain is overall spherical. Under different conditions, however, this picture may 
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break down.'^ For instance, if stretched by an external force or a net-charge repulsion, 
a PA chain will break into many smaller subunits. An important contribution to the 
PA free energy arises from the fact that PA charges on the periphery (the surface of 
the chain) are less effectively screened by other PA charges than those inside. This 
unfavorable free energy cost per unit area is called the surface tension, denoted as 
7. Following the scaling approach adopted in Ref.,*^ 7 is the number of blobs on the 
surface per unit area and is readily given by 



T 1 ifh f /b"^, low salt 



J2 



ifh) I {K-ob'^) , high sah 



see 



Finally, we have plotted the surface tension 7 as a function of kq in Fig. 4.6 
the inset). Since 7 is a direct result of PA effects, its magnitude is diminished as kq 
increases, as shown in the figure. 

While the notion of the surface tension will be useful in some contexts (e.g., the 
formation of a necklace globule!^, we only focus on overall spherical PA's here. An 
excess charge on a PA chain can drive a conformational transition to a pearl-necklace 
structure from a more spherical globule (referred to as a polyelectrolyte (PE) effect). 
This may be realized at a low-salt limit. It is worth noting that the low-salt limit 
here is different from that for the PA effect, since the PE effect is longer-ranged. 
Salt-screening is not felt sensitively by the PA effect, as long as kpa ~ fls/b'^ ^ «^o, 
which is A^-independent. On the other hand, the PE effect will be screened unless 
Rhq < 1 or < 1/-R. For the long-chain case, this imposes a very strict condition on 
salt concentrations. In conclusion, the low-salt limit for the PA effect can be realized 
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in a wide range of salt concentrations, while that for the PE effect is prohibitively 
narrow, as long as N ^ 1. We thus focus on overall "spherical" PA's, whether they 
are collapsed or swollen (e.g., self-avoiding walk chains). 

4.4 Polyampholytes with excess charges 

For several reasons, PA chains can carry an unbalanced, excess charge. Recall f± 
is the fraction of positive and negative charges on the PA chain. Then the chain 
will carry a net charge unless /+ = /_. The net charge repulsion now enters into 
the picture and competes with the PA effect. To focus on this polyelectrolyte effect, 
smear out all the charges on a PA chain, which is assumed to be inside an imaginary 
volume V. The excess electrostatic energy stored in this volume is expressed as 



Here the subscript of kq was dropped for simplicity; k is the Fourier conjugate to 
r — r'. Only k = contributes to Ep^, since the monomer density is assumed to be 
uniform inside V. The PE free energy thus scales as 





(4.27) 




B 



(4.28) 
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This estimate is relevant for the case R 3> It is tempting to relate this to the 



two-body term in Eq. 4.16 and interpret the PE term as renormalizing the second 
virial coefficient v, as in Ref.^ Strictly speaking, this reasoning is valid only if each 
"chain segment" remains spherical, in the sense that its width is equal to or compa- 
rable with its length. A crucial concept that describes the PE effect on chain shape 
is the electrostatic persistence length, denoted here as le, which was originally in- 
troduced over three decades ago.'^ While this concept has been widely used in the 
literature, ironically, there has been a controversy over whether Ig ~ /t"^, as originally 
predicted'^'^^'^ (also see references therein). 



With this subtlety in mind, let's insist on using the PE term in Eq. 4.28 and 



compare it with the PA term in Eq. 4.15 These two effects are comparable if FpE 



|Fpa|. For the case kq > kpa, this implies that 

+ (4.29) 

(One can arrive at this by considering Eq. (4.11) in Ref.l^. For a random PA chain, 
/+ = /_ and I/+ — /-I ~ 1/\/N. Despite some uncertainty in the k dependence of 



Eq. |4.29[ it is clear that the PE effect can be easily dominated by the PA effect. For 
different solvent chemistry or polymerization processes (i.e., /+ ^ /_), however, the 
PE term can be dominant and will expand the chain. As a result, R ~ A^'^, where 
z/ 3/5 if Rk ^ 1 or = 1 if Rk < 1. We will not attempt to refine earlier results 
for PA chains with excess charges, especially in regard to the k dependence of their 
size. As more results become available, they can be incorporated into our approach. 
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4.5 Elastic response of flexible polyelectrolytes 



The response of flexible polyelectrolytes such as ssDNA and RNA to mechanical force 
has further given a fundamental understanding of elasticity of PEs and has further 
clarified the changes in persistence length as the ionic concentration is varied. In fact, 
double stranded DNA molecules, which we have seen are intrinsically stiff polyelec- 
trolytes, have long been the major polymers of interest in single molecule stretching 
experiments.'^^'^ For dsDNA and related stiff PEs, a wormlike chain (WLC) model 
readily explains the measured force-extension data'^^'^^ as well as the distribution 
of the end-to-end distance.^Ql For flexible polyelectrolytes, however, it has been diffi- 
cult to account for the data using polymer models alone, thus raising the possibility 
that microscopic structures might matter. For these systems, it is probably the 
strong interplay between the electrostatic effects and the small intrinsic stiffness of 
the molecules that complicates the physics of stretching. For ssDNA and RNA in 
particular, base pairing between pairs of nucleotides along the chain backbone could 
be an additional complicating factor. Here, we will just focus on the situations in 
which the base pairing interactions are negligible. 

We begin with a survey of experiments to indicate the diversity of responses of 
charged flexible PEs. Dessinges et alJ^ used magnetic tweezers to obtain force- 
extension curve (FEC) for a 11 kilo base ssDNA molecules, in the range of force from 
O.OSpN to about 50pN. In order to extract structural information of the molecules, the 
authors used an extensible freely-jointed chain with electrostatic interactions. The 
electrostatics is modeled using the Debye-Hiickel potential with an effective charge 
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density, u, along the backbone first described by Zhang et alJ^ 



Eelec. = — I dSidSj ^ (4.30) 



where e = 80 is the dielectric constant of water, is the spatial distance between 
two points i and j, and the double integral is taken along the chain contour. The 
force-extension data for this model is determined through Monte Carlo simulations. 
It appears that the model reproduces well stretching data in several different ionic 
conditions. In particular, in solution of lOmM phosphate buffer with A^i = 1.87nm, 
where the base pairing interactions are suppressed, the FEC of the ssDNA appears 
to be almost straight for at least two decades of force in the log-linear scales (see 
Fig. 4 of ref.^ At the values oi u = 1.28e/nm and Kuhn length of 1.6 nm, the 
theoretical FEC follows very the experimental data. The extensible wormlike chain 
and the extensible freely jointed chain models without electrostatic interactions can 
only fit the FEC for forces at least 20pN or even higher. Even so the persistence 
length obtained with the extensible wormlike chain model is too small, Ip = 0.21nm. 
(It should be noted that although the FEC portion for / > 20pN was plotted in the 
same figure, the data had actually been obtained previously for overstretched dsDNA 
in a solution at much higher ionic concentration of 150mM Na+ by Rief et al.^^) 

In 2004, Seol et alJ^ built synthetic RNA constructs made of only uracils, or 
poly(U) (no possibility for base-pair interactions) and examined their elastic proper- 
ties using optical tweezers. The solutions considered have the ionic strength ranging 
from 5mM to 500mM, and the force range is from about 0.5pN to 50pN. The authors 
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used, instead, a modified extensible WLC model with electrostatic interactions where 



where S is the stretch modulus in unit of force (pN) and the scale dependent persis- 
tence length is Ip = 10 + lospf^iq)- Thus, it is assumed a priori that the electrostatic 
persistence length is given by the OSF theory. The modified FJC model with electro- 
static interactions could not fit the data well at high ionic concentrations. However, 
WLC model appears to reproduce the experimental FEC at high concentration of 
Na"*" down to about lOmM. The resulting persistence length does not depend on the 
force scale for concentration up to 500mM, whereas it decreases by two fold from 
/ = IpN to about 50pN. 

More recently, Saleh et al.'^^'^used magnetic tweezers to stretch ssDNA molecules 
that are specially treated to avoid base pairing even at the lowest forces. The ionic 
concentrations were varied over a broad range (20mM to SOOOmM of Na^), and the 
minimum force is about O.OSpN, whereas the maximum force is similar to those in the 
aforementioned works. Instead of modeling electrostatic through the Debye-Hiickel 
potential, they use the blob picture and scaling argument s'^^'^ for self-avoiding chains. 
Besides, demonstrating for the first time the Pincus regime x ~ at salt concentra- 
tions up to 2000mM, the data also shows the logarithmic dependence of the extension 
on the force, as observed earlier .fl^IllIlSI xhe two regimes are demarcated by a crossover 
force fc at a characteristic extension Xc- When all the FECs at different salt concen- 



the relative extension is given by'^'^ 




(4.31) 
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trations are normalized by fc and Xc, they collapse fairly well onto a single master 
curve. The scaling arguments show that ~ ksT/l and Xc ~ (f//'^)l/3, where / is 
the Kuhn length, which is twice the total persistence length, and v is the excluded 
volume between the Kuhn segments. Since fc and Xc can be extracted from the FECs, 
V and / then can be estimated without any assumption on their relation. As a result, 
/ ~ j-o.5i±o.4^ which is almost close to /^^-^ or the Debye length A^. In addition, 
at / = 3000mM, there is effectively no excluded volume effect due to the counterion 
condensation (see next section) that renders the excluded volume v negligibly small 
(O— condition). Under this condition, the FEC can then be fitted very well using 
the standard WLC modell^ but not with the FJC model. The extracted persistence 
length Ip is 0.62 ± 0.01 nm, which can be treated as the "bare" persistence length of 
ssDNA molecules. The seemingly diverse responses of a single chain can be explained 
using a unified theory (see below). 

In 2006, Toan et al.'^ used their Thick Chain model, which views any polymer 
as a tube with uniform thickness, to fit the FEC data in ref.'^ The excluded volume 
effects along the chain backbone leads to an effective persistence length. The extracted 
parameters are the effective monomer length and the polymer thickness, both of which 
can be combined in a simple formula to produce the persistence length. Besides 
producing a good fit to the data for almost all salt conditions, the model shows 
that the effective thickness of RNA varies linearly with the Debye length. More 
importantly, it has been shown the resulting persistence length can be fitted with 
/p = /O + cA/), where ZO ~ 0.66nm and c = 0.43 is a constant. Thus, in contrast to the 
results by Yeol et al.,^^ which used the WLC model with electrostatic interactions that 
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assumed a OSF type persistence length, the dependence here is hnear or somewhat 
sub-hnear. 

In a recent paper pubhshed in 2010, Toan and Thirumalai'^ using a combination 
of geometrical arguments and scaling theory to derive a model-free unified theory 
for semiflexible polymer stretching at high force. The theory can be used for poly- 
electrolytes in moderate to high salt concentrations. One of the main results of the 
theory is that the apparent elasticity of a polymer is inherently force-dependent. The 
FEC will appear to be that of a WLC model in the force range of = ksT/lp to 

~ AkBTlp/b"^ with b being the monomer length, and the polymer behaves as a 
FJC for forces higher than J^h- Although similar results for WLC-like models were 
obtained previously ,^2'^^ the results are general and hold good for any semiflexible 
chain regardless of the details of the monomer-monomer interactions. The numer- 
ical values of the two crossover forces for typical ssDNA and RNA molecules are 
J^h ~ 4pN and J^h ~ 50pN, which turn out to be relevant to the reported discrepancy 
in the elastic behaviors of structurally similar molecules, ssDNA and RNA. The higher 
crossover force is about the maximum force than can be achieved in most of the ex- 
periments,'^^'^^'^ which means that the WLC behavior is more likely to be observed. 
This is also true for the data cited in ref.,'^ where the FJC fits well the data above 
50pN (in fact the data is from another workl^ done at different ionic concentration 
(see above). In fact, fitting the same high force data using a generalized formula for 
both the WLC and FJC regime, Toan et al.'^show that the fit is reasonably good and 
produce a physically meaningful persistence length of 0.72nm. Moreover, the theory 
naturally leads to a simple, called the 1/2-rule, which states that the extension x at 
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force f = J-'i = ksT/lp is 1/2. That means the persistence length could be quickly 
estimated from the FEC without doing any fit, by calculating the force at x = 1/2. 
Indeed, the rule does lead to the estimate of the persistence lengths of RNAf^ and 
ssDNAfl^ in different ionic solutions, that are similar to each other and also scale 



linearly in the Debye length (see Fig. 4.5). 



4.6 Simulations of collapse of an isolated polyelec- 
trolyte chain 

A theoretical description of PE collapse is relevant in a number of applications includ- 
ing RNA folding and conformational fluctuations of disordered proteins. The initial 
event in the folding of RNA is counterion driven collapse. The kinetics of collapse of 
this process is driven by a number of factors such as solvent quality, valence, shape, 
and size of counterions. Although the flnal structure is unaltered the pathways in 
the coil-globule transition is dependent on the details of interactions between ions 
and RNA. In addition, DNA collapse into toroidal structures can also mediated by 
counterions. These considerations make the study of collapse of PEs (and the related 
PAs) important. 

Lee and Thirumalai^Sl analyzed the collapse of flexible polyelectrolytes in poor 
solvents using simulations. In their model, each monomer is approximated as a van der 
Waals sphere with radius b/2 that carries a charge of — e. Two successive monomers 
are connected by a spring with spring constant SksT/b"^. The polyelectrolyte chain is 
in a box with counter ions, whose valence varies from +1 to +4. The number of the 
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ions is such that the total system is neutraL The non-electrostatic interaction between 
the particles (monomers or counter-ions) i and j with radii (Xj and aj respectively that 
are at a spatial distance Vij is modeled using Lennard- Jones potential 



12 



(4.32) 



where tq = (crj + crj)/2. The parameter eLj is used to control the quality of the solvent, 
which is expressed in terms of the second virial coefficient V2 roughly equal to the 
excluded volume 



V2 



dh (1 



(4.33) 



For poor solvents or hydrophobic V2 < 0, and the strength of the hydrophobicity 



depends on the value of V2 (see Fig. 4.7). The number of monomers considered 
was = 240, for which the effective G-temperature corresponds to eLj = 0.5 and 
V2 = —1.96'^. The Coulomb potential is used to model interactions between the 
charged particles i and j with valences Zi and Zj respectively 



He {rjj) _ hziZj 
ksT r 



(4.34) 



The chain dynamics is followed using Brownian Dynamics simulations. The time unit 
is r = b'^(/2kBT, where ( is the friction coefficient of a monomer. 
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4.6.1 Effect of Valence and Solvent Quality on Collapse Dy- 
namics 



The dynamics of collapse of the PE chain is determined by a balance between hy- 
drophobic interactions and charge renormalized electrostatic potentials. The chain 
would be in extended conformation if the electrostatic repulsion is dominant. On 
the other hand, the valence-dependence effective electrostatic attractions between 
monomers due to the Manning condensation (of counterions onto the chain) could 
make the chain collapse. In particular, let us consider the case of near and hy- 
drophobic conditions. 

In near B solvent (eLj = 0.3 and V2 = —0.066^), the collapse kinetics is monitored 
through the normalized radius of gyration Rg{t)/Rc, where Rc — 4.086 is the size 
of compact globule. The collapse is strongly dependent on z, and occurs readily for 
z = 3 and 4 within the time of 150r from initial condition. For z = 2, the collapse 
only happens after 600t, whereas it is not seen for ^ = 1 even up to the time limit 
of the simulation. There are two effects of counterion condensation: (i) the overall 
charge of the polyanion is greatly reduced and (ii) there is an effective attractions due 
to the excess charges which effectively makes the solvent poor and gives rise to the 
collapse. In the hydrophobic condition V2 < —10.36^, on the other hand, the neutral 
chain would adopt a compact conformation. Thus the counterion condensation would 
further accelerate the collapse. In deed, with eLj — 2.0 and V2 — —25.86^, the collapse 
already occurs a,t t — 200t ior z — 2 and at about 600t ior z — 1. 

The collapse is also found to happen in three stages. The first stage corresponds to 
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the condensation of the counterions to the charged backbone, that is driven mostly by 
diffusion and occurs on the order of 25t independently of the valance z. The second 
stage is the formation of pearl- necklace structures of globular clusters containing both 
monomers and counterions. The clusters are mainly local, namely the monomers in 
a cluster are predominantly neighbors, and are connected by strings. In the third 
stage, all the clusters merge and the largest cluster grows at the expense of smaller 
ones (see below for a theoretical explanation) . 



The collapse mechanisms depend on the valency of the counterions. For z = 1, the 
monovalent counterions when condensed combine with the charges on the backbone 
to form random dipoles of magnitude p ~ eb. When the attraction between two 
dipole exceeds /c^T, i.e. the spatial distance is less than Ib, contacts between distinct 
segments of chain can form. Because this is a very short distance, the globular clusters 
are predominantly local and their sizes are very small while the number of them is 
large. Thus it takes very long time to reach the global compact structure. When 
the counterions are multivalent, they can both neutralize the negative charge on the 
monomer and provide excess charge of \{z — l)e|. Bare monomer charges that are 
separated by a large distance along the contour can be attracted to the positive charge 
in a process called "ion-bridging" . The range over which such attractive interactions 
are effective increases with z. Thus the size of the initial clusters and the efficiency 
of the collapse also increase with z. 
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4.6.2 Phase diagram in the {v2, z) plane 

The observations concerning the structure of the collapsed globules together with 
simple scaling arguments can be used to construct a valence-dependent phase diagram 
for strongly charged polyelectrolytes in poor solvents {v2 < 0). When counterion 
condensation takes place, the decrease in the effective charge of the polyanion can be 
computed by equating the chemical potential of the free and condensed counterions 
(two-phase approximation). Generalizing the arguments of Schiessel and Pincus'^to 
arbitrary z, the authors find that the total charge of the polyelectrolyte decreases 
from Nfe (/e is the fraction of charge, which is 1 here) to Nfe ~ k{L/lB){l/z) where 
k ^ — In and (p is the volume fraction of the free counterions. The size of the 
polyelectrolyte is given by L ^ k'^lP'N /Igz'^ provided Ib > k'^z~'^bN^^'^. As the quality 
of the solvent decreases to a level such that thermal blob size ~ b'^/\v2\ < C,ei ~ 
iBz'^/k'^ (size of the electrostatic blob) then the chain condenses to a globule. The 
boundary dividing the stretched and collapsed conformation is obtained by equating 
^ei and S,T and is given by \v2\ ~ b^k'^/lBz'^. 

In the globular phase, two regimes are found, one corresponding to the Wigner 
crystal and the other a Wigner glass. The boundary between the two is obtained 
by equating the gain in the energy upon condensation ze'^/de) to the attractive 
interaction due to the poor solvent quality {kTv'^b~^). This leads to the condition 
\v2\ ~ z^^'^b^^'^l^J'^ . To determine the boundary more precisely we have to account 
for the induced attraction between monomers after counterion condensation. Thus, 
the actual determination of the boundary separating the Wigner crystalline regime 
and the glassy regime requires equating ze'^/de and fcTf|^6~^ where V2r is the renor- 
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malized second virial coefficient, and d is approximately the distance between the 
condensed counterion and the backbone charge. This argument shows that the shape 
of the PE chain requires, in a nontrivial way, the couphng of electrostatic interactions 
and effects coming from solvent quality. In light of this, the boundary indicated in 



Fig. |4.7| should be regarded as qualitative. 

The boundary between the two collapsed phases is difficult to determine from the 
simulation because of the extreme slow dynamics in the glassy phase. For the purpose 
of illustrating the phase diagram, we assume that the ionic glass is an equilibrium 
phase. To validate the phase diagram, extensive simulations were performed for all z 



1, 2, 3, 4) at values of 1^^ and \v2\ indicated by asterisks in Fig. 4.9 The structural 



conformations are in qualitative accord with the predicted phase diagram. 



4.7 Theory of Collapse Dynamics 

Although there are a number of theoretical models that describe the collapse kinetics 
of homopolymers with applications to protein folding there is only one study that 
has considered the rate of PA and PE collapse. The theory proposed by Lee and 
Thirumalai'^ (LT) adopts an approach developed by Pitard and Orlandl^ to describe 
uncharged homopolymer collapse in poor solvents. For simplicity, we consider the 
collapse of PAs, which as will argue bears a strong resemblance to the behavior 
expected for PEs as well. 

A model system, which may captures some generic aspects of biomolecules, is 
a polyampholyte, which is a linear polymer chain that has both positive and neg- 
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ative charges along the backbone. The very large muscle protein titin associated 
with sarcomere contains regions that have PA characteristics as do many intrinsically 
disordered proteins. Because of the simultaneous presence of positive and negative 
charges the shape of PAs is determined by competition between the two conflicting 
interactions. The repulsion between the charges tends to swell the chain (the poly- 
electrolyte effect) while the attraction (polyampholyte effect) tends to collapse the 
chain. It is known that if the total charge, Q, of PA is less than ~ ±Ne {N is 
the number of monomers) the chain adopts a globular conformation, and is extended 
otherwise. If Q > Qc then it can be shown, using an analogy to Rayleigh instability 
of charged droplets, that the chain can be visualized in terms of pearl-necklace struc- 
tures, which is similar to that found in weakly charged PEs. The LT theory showed 
that collapse to a global globular structure, which can only take place occur when 
Q < Qc, occurs in two major steps. In the first stage the chain forms a metastable 
pearl-necklace structure reminiscent of the equilibrium structures for Q > Qc- In the 
second stage the pearls (domains) merge leading to the compact globular conforma- 
tion. The theory for the intermediate time regime (when pearl-necklace structures 
form) was developed by adopting the procedure suggested by Pitard and Orland.'^ 
The LT theory showed that the square of the radius of gyration decays as 



Rlit) ^ Rimi - {—)r (4.35) 

where a = | and the characteristic time tpa ~ A^^. On the time scale tpa metastable 
pearl-necklace structures (local structures connected by strings) form. Interestingly, 
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formation of such structures have also been observed in the coil-globule transition in 
uncharged homopolymers. 

Dynamics of PE collapse: The theory described for PA can qualitatively predict 
the collapse dynamics of strongly charged PEs. Let a strongly charged PE (fraction 
of charged monomers is close to unity) chain be initially in the 6- solvent with respect 
to the neutral polymer. At t — 0, we imagine a quench to low enough temperatures 
so that counterions start to condense. Upon condensation of counterions the confor- 
mation of the chain approximately resembles that of a PA. The relevant length scales 
for PE are the electrostatic blob length D Ki and ~ clq^^ where z is the 
valence of counterions, k = —Incj) {(p is the volume fraction of free counterions), and 
6 is the collapse temperature of uncharged polymer. If is not too small, the PE 
chain undergoes a sequence of structural changes en route to the collapsed conforma- 
tion. After counterion condensation PE evolves toward a metastable pearl-necklace 
structure. For weakly charged polyelectrolytes such structures are the equilibrium 
conformations in poor solvents. The dynamics of this process can be described using 
the theory developed to describe collapse of PA. We assume that shortly following 
the quench to low temperature the counterions condense onto the PE chain. The 
time scale for this process is diffusion limited, with values that are much smaller than 
time in which the macromolecule relaxes. Upon condensation of a multivalent cation 
{z > 2), the effective charge around the monomer becomes {{z — l)e). If the locations 
of the divalent cation are random and if the correlations between the counterions are 
negligible then in the early stages, a PE chain with condensed counterions may be 
mapped onto an evolving random PA. With this analogy we suggest that the pearl- 
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necklace-like structures that are found in the collapse of strongly charged PEs should 
also be governed by Eq. (35). 

Late stages of collapse: At long times the pearl-necklace structures merge to 
form compact collapsed structures. This occurs by the largest cluster growing at the 
expense of smaller ones, which is reminiscent of the Lifshitz-Slyozov mechanism.'^Zl If 
this analogy is correct then we expect the size of the largest cluster S{t) to grow as 
S{t) ~ t^. The collapse is complete when S{t) ~ Rg{t — ?■ oo) ^ which implies 
that the characteristic collapse time is ~ iV. 

Langevin simulations for strongly charged flexible PEs where the collapse is in- 
duced by multivalent counterions have shown that late stage coarsening indeed occurs 
by the Lifshitz-Slyozov growth mechanism. Figure 1 shows a representative example 
of the conformations that are sampled in the dynamics of approach to the globular 
state under B-solvent conditions. Both panels show that in the later stages of collapse 
the largest clusters in the necklace-globule grow and the smaller ones evaporate. This 
lends support to the proposed Lifshitz- Slyozov^^ZI mechanism. To estimate , we have 
used simulations to calculate the number of particles that belong to the largest cluster 
Ns{t) as a function of time. Figure 2 shows the linear increase of of the number of 
particles Ns{t) at long times from which the growth time increases linearly with N 
as expected from the LS mechanism. . The change of slope for long times is due to 
the finite size effects and indicates the completion of the globule formation. 

The mechanism of approach to the globular state for PA and PE should be simi- 
lar. The collapsed conformations are reached via metastable pearl-necklace structures. 
For PE the driving force for forming such structures is the counterion mediated at- 
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tractions. Charge fluctuations in PA lead to pearl-necklace structures. The lifetime 
of pearl-necklace structures depends on energy barriers separating merged and un- 
merged clusters. The necklace-globule conformation consists of n globules with nearly 
vanishing net charge that are in local equilibrium. The free energy of the i^'^ glob- 
ule is Fi ~ (^)AF_R? + A-naR} where Ri is the radius of the i^^ globule and a is 
the surface tension. Note that A/ is the same before and after the merger of clus- 
ters. The free energy difference between a conformation consisting of two clusters 
and one in which they are merged is AF ~ (Svrcr — 4(2)37rcr)(^)i. Charge fluctua- 
tion in each globule is qo{N/ny^'^. If the Coulomb energy fluctuation of each globule 
{6E ~ [N /n)q^/ aQ^^)'^) free energy difference between the conformation with two 
separated clusters and the one where they are merged, then the system spontaneously 
grows to a large cluster. 

4.8 Conclusions 

We have only described a few topics dealing with properties of isolated polyelec- 
trolytes and polyampholytes with some applications to biological systems. There are 
a variety of issues especially such as the viscosity of PEs that are puzzling and war- 
rant explanation. In addition, bundling of charged systems (actin and microtubule for 
example) is another topic with applications in biology that have not been adequately 
described. Only recently, polyelectrolyte concepts have been applied to understand 
conformational fluctuations in intrinsically disordered proteins.'^ Detailed theoretical 
and experimental studies are needed to clarify the role that polyelectrolyte effects and 
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Figure captions 



Figure 



4.1 



(Reproduced from'^ The dependance of R^ on Mg^"*" (squares) and 
Na"*" (triangles) concentration. Solid lines are the fits using Hill equation in the form 
A{1 — [Mg^^]"/ (C^ + [Mg^"*"]") } + ?/05 where A, n, and yo are adjustable parameters. 
The best fit values for n are 3.33 and 1.20 for Mg^+ and Na+, respectively. Clearly, 
cooperativity of the folding transition increases as the valence of ions increases. 

(Reproduced from^^ Distance distribution functions P(r) at various 



Figure 



4.2 



Na"*" concentrations at 32 °C are obtained by inverting the measured scattering inten- 
sity by a Fourier transform. The solid lines are fits of P(r) ~ exp {—jz^), which is 
the expected asymptotic result for WLC. The concentrations of counterions are given 
in the insets. 

Figure 4.3: (Reproduced fiou^^ Dependence of Ip on 1/k^ in Mg^"*" (solid circles) 
and Na"*" (solid squares). Lines represent fits to the data. Note that the 1/k^ for Mg^+ 
is given on top. 

Figure 4.4 (Reproduced fromf^l) Variations of the experimental persistence length 
with the reciprocal of the square root of the ionic strength for sodium poly [( (aery lami- 
do)methyl) propanesulfonate] (▼), sodium (carboxymethyl) cellulose (□), sodium poly- 
acrylate (o), sodium poly(styrenesulfonate) (•), and the sodium salt of an alternating 
copolymer of isobutyl vinyl ether and maleic anhydride at two degrees of neutraliza- 
tion a = 0.5 (A) and at a = 1.0 (V). 

Figure 4.5 (Reproduced from^*^) Analysis of the experimental results. Estimates 



of Ip{Xd) for poly(U) (circles) and ssDNA (filled symbols) using the 1/2-rule. Linear 
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fits to the data (solid and dashed fines) yielded tfie bare persistence lengtfis Ip ~ 
0.67nm for poly(U) and Ip ^ 0.63nm for ssDNA. 

Figure |4.6 The equilibrium size of a polyampfiolyte cfiain as a function of 



for a few cfioices of /: JIb = 0.1,0.2,0.3. Our results sfiow fiow tfie PA cfiain 
is swollen by salt-screening, wfiicfi tends to weaken tfie polyampfiolyte effect: tfie 
competition for tfie screening of PA cfiarges by salt ions is more pronounced at fiigfi 
salt concentrations. Also plotted is tfie size of polyampfiolyte blobs, ^pa (tfie same 
color scfieme used for R-kq grapfis is adopted). Our result for R and ^pa are useful for 
illustrating tfie crossover from tfie low-salt to fiigfi-salt limit. Finally, we fiave plotted 
tfie surface tension of a PA cfiain, 7 ~ I/'Cpa ^ function of (see tfie relevant 



discussion around Eq. 4.26) 



Figure 4.7 Second virial coefficient as a function of eLj from Eq. (4.33). Tfie 



classification of solvent quality based on tfie values of V2 are sfiown. 



Figure 4.8: Tfie time dependence of tfie radius of gyration following a quencfi 
from G-solvent to poor solvent condition. Tfie value of = 240 and Ib = 5.36. Top 
panel is for V2 = —0.066^ wfiereas bottom panel is for V2 = —3.696^. Tfie numbers on 
tfie curves denote tfie valence of tfie counterions. 



Figure 4.9: Valence dependent diagram of states in tfie {\v2\ and /^^) plane for 
strongly cfiarged PE. Tfie dasfied lines represent tfie boundary between stretcfied and 
cofiapsed states and depend on z. Tfie z-dependent solid lines in tfie collapsed region 
separate Wigner crystalline region from tfie Wigner "glassy" region. Tfie asterisks 
are tfie simulation results witfi eacfi data point corresponding to four z values {z = 1 
- 4). Pictures of tfie conformation of tfie cfiain in region A,B and C are also sfiown. 
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The l^^ values are indicated by arrows. The letter (a-1) near the asterisks corresponds 
to values of \v2\/b^^ 62.34,62.34,25.82,25.82,15.20,12.36,7.51,4.47,3.91,3.69,0.34,0.06 
respectively. 
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